Physical solution of ID quantum A^-boson system with zero-range pair interaction 
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The mathematically exact solution of a one-dimensional (ID) quantum A'^-identical-boson system 
with zero-range pair interaction has been well known. We find that this solution is non-physical, 
since there exists a paradox of its energy expectation value, leading a basic contradiction to the 
mean-field theory of the system. A new integral equation that is equivalent to the corresponding 
Schrodinger one is established and its formally physical solution is derived. Energy correction from 
the average potential and harmonic waves of different momentums are demonstrated, and scattering 
amplitude of the physical solution and mean-field theory of the system are discussed. 



m : 
o . 
o ■ 

<N ■ 

Oh' 
(D ■ 

00 : 



PACS numbers: 03.65.Ge; 03.65.Nk; 05.30.Jp; ll.80.Fv 



On 

o 
m 
o 

a ■ 

a • 

^ : 



The quantum integrability is an essential property of 
quantum mechanics and field theory. A kind of non- 
trivial integrable Hamiltonian operators, which describes 
some systems consisting of N identical non-relativistic 
particles interacting pairwise through some forms of po- 
tentials has been generally demonstrated by using differ- 
ent methods The forms of potentials contain the 
zero-range infinite pair interaction and the long- 
range finite pair interactions Although the inte- 
grable models may not be exactly solvable, the exactly 
solvable models are included by them. Much interest 
have been focused on the exactly solvable models, be- 
cause of their applications to the Bose-Einstein condensa- 
tions quantum Yang-Baxter equation (^], quan- 
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|22|. quantum Hall liquids (23j , , .24] . soliton theory |25| . 
and the black holes |2^ . 

The infinite zero-range pair interactions can be ex- 
pressed as the S function potentials, which approximate 
some practical short-range interactions and seem to be 
easier for mathematical treatment than the latter. How- 
ever, singularity of the d function may bring us some 
new troubles that lead to non-physical results. This can 
be shown through the simplest example, the Schrodinger 
equation Hip = —■^ipx^xr + c5(xr)ip = Eip in atomic unit 
and with the 5 function potential c5{xr)j where Xr de- 
notes the relative coordinate of two particles, ip{xr) rep- 
resents the relative wave function, c > is the zero-range 
interaction intensity and E the eigenenergy of the Hamil- 
tonian H . The standard quantum mechanical book 27\ 
had mathematically given its general solution as the lin- 
ear combin atio n of plane waves, ip = ae^'^^'' + l)e~^'^'^'' 
with k = \J1E and a, h being constants determined by 
the normalization and boundary conditions. This result 
had also been extrapolated to A^-body zero-range pair 
interaction problems and been recently applied 

to few-particle Bose-Einstein condensates |2^. However, 
there is a paradox of the energy expectation value in 
this mathematical solution. Noticing that the function 



giifcx^ denotes the state of a free particle, in this state 
the average kinetic energy is equal to the total energy 
and the average potential is {ii\cb{xr)W) = c|V'(0)p. The 
nonzero average potential leads to the paradox of average 
energy, namely the energy value may be (^/^liJ = E or 
(tplHlip) — E + c\ip{0)\'^. Both results are in agreement 
only for the identical- fermion system with ip{0) — 0. 
While for an A^- identical-boson system with tp{0) 7^ 
the mean-filed theory of second quantization is based on 
the average Hamiltonian and therefore there exists basic 
contradiction from the above result. Although adopting 
the finite short-range interaction 14], 29J- 31] instead of 
the infinite zero-range interaction can avoid the physical 
confusion, the Gross-Pitaevskii's nonlinear Schrodinger 
equation comes from the application of the zero-range 
interaction, which has been successfully used to treat 
the Bose-Einstein condensates ^13J- 15j. This means ex- 
istence of the physical solution of ID A^-boson system 
with zero-range pair interaction. 

In the previous works [s^], H^, we estabhshed and 
employed an integral equation to find the complete infor- 
mation implied in Schrodinger's quantum systems. Ap- 
plying this idea to write down an new integral equation 
of the zero-range pair interacting system and seeking its 
formally physical solution to eliminate the contradiction 
on energy expectation value are our main motive in this 
paper. In the physical states, the energy correction from 
the pair interaction and harmonic waves with different 
momentums are investigated, and the scattering ampli- 
tude and self-consistent field theory of the system are 
discussed. 

We commence with A'^ slow identical bosons propagat- 
ing in x direction and assume that the pairwise interac- 
tion consists only of a hard core of ID diameter. For 
small diameter the system is conveniently treated as im- 
penetrable point particles with the delta- function pseu- 
dopotential cS{xj — Xi), where c > is the interaction 
strength, Xi, Xj for i, j = 1, 2, • • • , A^ are the coordinates 
of particles i and j. The stationary state Schrodinger 
equation pertinent to the system reads as Q-Q 
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with m and E being the mass of a single particle and the 
energy of the whole system respectively. Save for the col- 
lisions, which occur only for Xi = Xj, the previous works 
treated the system as N free particles and constructed 
its exact plane wave solution mathematically ^J-^J- Ob- 
viously, the above-mentioned paradox of energy expec- 
tation value exists in this solution. In order to seek the 
physical solution without the paradox, we must count in 
energy correction from the pair interaction and solve this 
equation by using a new approach. 

Under assumption of the pairwise interaction and im- 
penetrability of particles, the exchange symmetric wave 
function can be expressed as Q, 

p odd i 

= X! n "^(^P^ + l' ^P^)^iXp^+l -Xp^)■ (2) 

p odd i 

Here denotes a permutation of xi • • • a;Ar, ^ 

represents sum of all the permutations, 0{xp.^-^ — Xp. ) 
is the step function satisfying 6[xp^^-^ — Xp^ ) = 1 for 

Pi-\- 1. "^Pi ^'Hd ^("^Pi-j-l *^Pi ) ' — " ^ for Pi-^ X "^Pi ' 

and 4'{xp^^-^^, Xp.) is the wave function of two particles. 
To form whole particle pairs, we take the large number 
N as an even number and approximate |l3j|-|l5l| iV ± 1 
to the even N. The exchange symmetry of the identical 
bosons requires (/)(a:p; , a^p^+i) = (/'(xp^+i, Xp.). Substitut- 
ing Eq. (2) into Eq. (1), we get the Schrodinger equation 
of two-boson wave function [j] , 

- H ^ + c(5(xp,^, -Xp^)^^ E.cj), (3) 
j=i Pj 

where Ei is the energy of two particles, obeying Ei — 
E for i = 1, 3, 5, • • • , iV — 1. Introducing the mass-center 
and relative coordinates 

Xi(^ = —j={Xp^j^-^ -\- Xp^), Xir — "^^(^pi+i •^Pi ) (4) 

to Eq. (3) and setting 

4>{Xpi^^, XpJ = (t)f^{Xic)(l)fJ{Xir), Ei = Eic + Eir, (5) 



we obtain the equation of relative motion 

and the mass-center motion equation. In the calculation, 
the formula S{xp.^-^ — Xp.) = S{\/2xir) = 5{xir)/\/2 has 
been employed. The mass-center equation is that of a 

free particle with eiginenergy Eic and eiginstate cf)"^^ = 
^^gi/cic^ic _|_ ^^^-iki^xi^ ^ where a nd bi are undetermined 

constants and kic = \p2mEiiJW . 

The corrected energy from the delta potential is of the 
order c|(/)jy^ (0)p, which is set as E[^, so we have 

E„ = Ef^ + E'„, kl = 2mE„/h^ = k^f'^ + k'^^, 
(0) _ , _ l^mEi V2mc 

The corresponding momentum reads as hkij. « f>\kf} + 

k'^/(2k^^)]. Inserting them into Eg. (6) yields a simpli- 
fied equation, then following Ref. leads to the new 
integral equation 




X {e^'^^"'-- j e-^^'^'--[c'5{x,r)-kfM'':^dx„ 

- e-^'^.v^- j e'^'S}^^,-[c'5{x,r)-kZ^t^dx„], 

(8) 

where B, C and D are undetermined constants obey- 
ing C < 0, D < 0. Applying Eq. (8) to Eq. (6), we can 
directly prove the agreement between them. Particularly, 
the relative wave function (/),[^'' in Eq. (8) does not satisfy 
the free particle equation such that the above paradox of 
energy expectation value is avoided. Completing the in- 
tegrations in Eq. (8), we get its formally exact solution 
as 



A- 



c'0l:?^(O) , ^ki 



+0(-x,r){[A+^ I e-'=^^"■4^)dx..]e^'=^^--H[B-^ / 



B- 



2^k^ 2fc(°) Jd 



2fcL°) Jc 



2k Jd 



"■}, (9) 



where (j)^^' {0) is the value of (pi^' at Xir — 0. Applying Eq. (9) and the above-mentioned (pi^' to Eq. (5), then to 
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Eq. (2), we obtain the formally exact solutfon of the A''- 
identical-boson system. The formal solutfon, of course, 

(%>) 

is still an integral equation, since 0-^ appears in the in- 
tegrations of Eq. (9). It is worth noting that the integral 

equation (9) completely describes the continuity of (Z)-^"* 

and discontinuity of dij)-^ /dxir caused by the singularity 
of delta potential. 

When the constants fc^^, A and B are taken as — 0, 

A^B + ^<^rX^l B^A-^<f>^iO) (10) 

with kjj. = /c-°\ Eq. (9) becomes the well-known solu- 
tions 3, 



which implies 



0{Xir) 



(11) 



{0) = A + B. 



The selection to A and B in Eq. (10) is because of the 
requirement of exchange symmetric Eq. (11), (l)\^\xir) = 



/,(p) 



{-x^r)- Solving Eq. (10) and Eq. (12) yields 



Applying this to Eq. (11), the normalization of (/)^^^ gives 
(f>[^^ (0) oc 1/ Vl with L being the length of motion region 



adjusted by the boundary conditions f^. Obviously, for 
the plane wave solution (11) the average kinetic energy 
equates the total relative energy so that the nonzero av- 
erage potential leads to the paradox on energy expecta- 
tion value. Only for an infinite L the average potential 
c|(/)(0)p vanishes and the contradiction does not exist. 
It is the non-zero energy correction E'^^. that helps us to 
overcome the difficulty on energy expectation value for 
the case of finite motion region. 



When fc-^ ^ is set, Eq. (9) cannot be exactly solved. 
However, it allows existence of the infinite series solu- 
tion in terms of plane waves of different wave vectors. 
The series expansion leads the terms proportional to k'^^ 
to contain the harmonic waves of various momentums, 
which are generated by the pair interactions. 

For fcf"' > k'^^ we can treat the terms proportional to 
k'^ as perturbation and construct approximate solution 
of the equation easily, through the unperturbed solution 



(12) (11). This requires c\4r\o)\^ ~ K < since E^^ 



is the corrected energy from the average potential. Note 
that only for small c or big L (small (f>'^j^\o)), it is valid, 
because of the small E'^^. In fact, simply inserting Eq. 

(11) into left hand side of Eq. (9) to replace (j)^^^ , we 
immediately arrive at the perturbed solution [32|, [33 ■ 



Noticing ktr « A: -"^ -|- / (2kl''^' ) , for simplicity we let 
kir in denominator of Eq. (13) be equal to k^^^ approxi- 
mately. To avoid resonance, this approximation does not 
be used for the exponential functions in Eq. (11). The 
perturbed solution then reads as 



^(0)^ 



J 



Be' 



1 + e~'''i?'"*'-/(^''i^') - e~*''"-^/(^''»°'^ 



Ae 



Ae' 



-ik'ilxir/(2k\"') 



-ikZD/{2kfJ) 



Be 



-ik): 



-}+0{k: 



'2) 

ir J ' 



where A and B obey Eq. (13) with kir = kf^ , 0{k'^r) 
denotes the terms proportional to fc^^. The exchange 

symmetry of </)-^^ infers e*'^^^'^/^^'^'^'^ — e^*'^^^-^/^^'^*^'^ 
namely -/c,'2(C + i:i)/(2fc^^°^) = 2n7r for n = 1, 2, • • • . The 
terms containing /c^y give perturbed correction to unper- 
turbed solution (11) that includes two different waves 
with wave vectors k^^^ and kir respectively. 

To simplify the discussions, we rewrite Eq. (14) as 



(0) 



^i^^lfe' =0, and let Hi, 



be the Hamiltonian associated with Eq. (6). Thus the 
constant 4>f^ (0) is determined by the normalization of 



, while the relation among k[^, k'f) and c can be 



give by 



Hir\(t>. 



ip)\ 



e\^^ -t- E'jy, no paradox of the 



energy expectation value occurs. Therefore, the solution 
(14) is physically correct one. Noticing {(pfj \Hir\(j)"^) — 
£"1"' + c|(/)-°^(0)P/a/2, we have the corrected energy 

h^k'I 



EL 



^-^|4-^(O)|2 + (0^^|iJ„|4-^) 
{^fJ\H,r\^l) + {cj,l\H„\^l). 

(0)/ 



(15) 



Clearly, only for infinite L, 0-^ (0) — 0, Eq. (15) allows 



E^ = 0, k^ = such that 



a(p) _ a(o) 



, Eq. (11) becomes 
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exact solution of the integral equation (9). 

Given the physical solution (14), the scattering ampli- 
tude can be approximately obtained. We apply k'^^ « 
and Eq. (13) to Eq. (14), resulting in the scattering 
amplitude 

, \ _ fcir + ic' ^ k\^^ + k'-^/{2k\J) + ic' 

^''^ - B-k„- ^c' - 40) ^ ^./(24°)) - ^c' ^ ^ 

which is similar to the known result 3-0- This sim- 
ilarity infers that we can directly extend this result of 
two-particle problem to iV-boson system. 

It is impossible, of course, to seek perturbed solution 
of the integral equation (9) for large average potential 
that necessitates the mean-field approximation. After 
eliminating the paradox on energy expectation value, 
the contradiction implied in mean-field theory of ID A'^- 
boson system will be avoided. In fact, adopting Eqs. (2) 
and (9), the average Hamiltonian of second quantization 
reads as {i>\H\il}) with H including the kinetic energy and 
internal pair interaction, even external potential. Dif- 
fering from the plane wave solution, the formally exact 
solution (9) does not identify the average kinetic energy 
with the total average energy and allows existence of the 
nonzero average potential thereby. Thus the well known 
mean-field theory is valid for the zero-range pair inter- 
acting system that means the wave function obeying a 
nonlinear Schrodinger equation 12]- 15]. 

In summary, we have investigated the ID quantum iV- 
boson system with infinite zero-range pair interaction. 
The paradox of energy expectation value was found for 
the plane wave solution previously reported in some ar- 



ticles. The integral equation (8) that is completely on 
an equality with the corresponding Schrodinger equation 
was established and was applied to construct the formally 
physical solution (9) without the paradox. The approx- 
imate solution and energy correction from the nonzero 
average potential were detailed, and the scattering ampli- 
tude, harmonic waves of various momentums and mean- 
field theory of the system were discussed simply. 

The method and results of this work are useful for 
treating any quantum system with the singular delta po- 
tential. For instance, in the case of a delta potential well 
with c < 0, the relative energy Eir is negative and the 
"wave vector" kir = i\kir\ is imaginary. Inserting these 
into Eqs. (8) and (9) and selecting suitable constants 
A, B, C and D to satisfy the boundedness conditions as 
in Ref. js^] , we can directly construct the bound state so- 
lution of the system with nonzero energy correction and 
without any paradox. It is straightforward to extrapolate 
the results to the ID system with combined potential of 
the pair interaction and harmonic trap, which has often 
been employed to treat the Bose-Einstein condensates. 
It is worth noting that in the bound states the integral 
equation (8) has no physical exact solution even for in- 
finite motion region that necessitates approximate and 
numerical solutions. The quantum many-body problem 
is a source of the quantum Yang-Baxter equation that 
is worth the further investigations. 
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